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FOREWORD 


In Part I of this report (Schwiderski, 1978a), a combined hydrodynamical-empirical method 
was introduced to compute numerically harmonic partial tides in the world oceans with an ac- 
curacy of better than 5 cm, which is needed in various military and civil applications of today. In 
this report, the computed semidiurnal principal lunar tide (M, ) is displayed in an atlas of tabulated 
tidal charts and plotted corange and cotidal maps. 
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ABSTRACT 


In Part I (Schwiderski, 1978a) of this paper, a unique hydrodynamical interpolation technique 
was introduced, extensively tested, and evaluated in order to compute partial global ocean tides in 
great detail and with a high degree of accuracy. This novel method has been applied to construct the 
leading semidiurnal principal lunar (M, ) ocean tide with an accuracy of better than 5 cm anywhere 
in the open oceans. The resulting tidal amplitudes and phases are tabulated on a 1° x 1° grid system 
in an atlas of 41° x 71° overlapping charts covering the whole oceanic globe. A corresponding atlas 
of global corange and cotidal maps is included to provide the reader with a quick general overview 
of the major tidal phenomena. The specifying hydrodynamical parameters of the model are listed 
along with quoted sources of empirical tide data, and significant tidal features are explained and 
discussed. 


1. INTRODUCTION 


Part I of this paper (Schwiderski, 1978a) introduced a unique combination of hydrodynamical 
and empirical methods to model detailed ocean tides with a component accuracy of better than 
5 cm anywhere in the open oceans. This enormous accuracy is well above minimum requirements 
set by, for instance, the National Aeronautics and Space Administration (NASA) and the Depart- 
ment of Defense (DoD) — to map the geoid at sea by satellite altimetry to within 10 cm. The fol- 
lowing features of this unique hydrodynamical interpolation model made the achievement of this 
accuracy possible. 


a. A spherically graded 1° x 1° grid system is set up in connection with a corresponding 1° x 
1° bathymetry to assure a sufficient resolution of all important tidal phenomena. 


b. The bathymetry of the gridwise, simply connected ocean basin is hydrodynamically defined 
(Schwiderski, 1978c) by appropriate modifications of earlier realistic depth data collections. The 
hydrodynamical redefinition was needed in order to model the well-known strong distortion and re- 
tardation effects of shallow continental shelves, narrow ocean ridges or island chains, and other sig- 
nificant bottom irregularities (see tidal charts in Appendix A). 


c. The Boussinesq substitution of the turbulent Reynolds stresses is applied in the form of 
eddy dissipation with a novel physically meaningful eddy viscosity that depends linearly on the 
lateral grid-cell area and, hence, directly on the ocean depth. 


d. The linear law of bottom friction is introduced with a bottom-friction coefficient depending 
linearly on the bottom grid-cell area which is independent of the ocean depth. In boundary cells, 
the otherwise constant friction coefficient is subjected to an indirect cellwise adjustment in order 
to permit a consistent hydrodynamical interpolation (see h., below) of empirical tide data known 
from tide gauge stations at continental shores, islands, or other shallow-ocean bottom irregularities. 


e. The effects of the terrestrial tide and the oceanic tidal load are included as simple second- 
order approximations in the sense of Love and Accad and Pekeris (1978). 


f. The Hansen-Zahel (Zahel, 1970 and 1977; Estes, 1977) finite differencing technique is 
modified by anew differencing scheme in time which improved decay, dispersion, and stability char- 
acteristics of the numerical procedure and facilitates the simple indirect adjustment of the bottom- 
friction coefficient in the hydrodynamical interpolation technique (see d. and h.). 


€ 
g. At land-ocean cell walls, the conditions of no-flow across and free-slip along the boundaries 
are enforced. The no-flow condition is subsequently relaxed by allowing controlled periodic inflows 
and outflows over the mathematically assumed boundaries. This allowance redefines indirectly more 
realistic shorelines in order to further improve the consistency of the hydrodynamical interpolation 
of empirical data (see d. and h.). 


h. A unique hydrodynamical interpolation technique is introduced which incorporates into the 
theoretical model empirical tidal constants collected from over 2 000 tide-gauge stations around the 
world in a hydrodynamically consistent fashion (see d., f., and g., above). 


soa With these features, the new model was successfully applied to chart the semidiurnal principal 

lunar (M2 ) ocean tide with the desired accuracy. The technique and accuracy of the model were ex- 

| tensively described and discussed in Part I of this paper as well as in subsequent journal publications 
and symposia presentations by the author (Schwiderski 1978a, b; 1979a, b, c,d; and 1980). 


{ During the past year, after the construction of the first M, tidal chart, the author received in 
private communications numerous new or improved empirical ocean-tide data from Gill and Porter 
(1978, from deep-sea stations in the New York Bight), Wyrtki (1978, from island stations in the 
Pacific Ocean), Zetler (1978, from deep-sea and shore stations in the southeastern Indian Ocean; 
see Irish and Snodgrass, 1972), and Cartwright (1979, from 108 deep-sea stations in the North 
Atlantic, northeastern Pacific, and southeastern Indian Oceans; see Cartwright et al., 1979). In ad- 
dition, in early 1979 the author received the entire (3 898 stations) worldwide ocean-tide data col- 
lection of the International Hydrographic Bureau (1978), which was updated and taped by the 
Canadian Marine Environmental Data Service. 


With few exceptions of apparently erratic tabulations, the empirical data from ali those sources 
substantiated everywhere the estimated enormous accuracy of the tide model. For example, the 
deep-sea profile measurements of tides in the New York Bight by Gill and Porter (1979) verified the 
globally modeled, incredibly sharp tide decay (Figures a and b) across the continental shelf of less 
than 200-km width from a range (= double amplitude) of 4.6 feet at the coast to an almost con- 
stant range of 3 feet in the deep ocean (see Table Sa in Part I (Schwiderski, 1978a) and the new 
slightly improved Tables 8M or 8N in Appendix A, which include the Gill and Porter (1978) data in 
interpolated form, indicated by subbrackets). 


For the novel hydrodynamical interpolation technique, it was emphasized in Schwiderski 
(1978a, 1979d, and 1980) that hydrodynamically interpolated marginal empirical tide data in iso- 
lated grid cells distort the surrounding computed ocean tide only insignificantly. This fortunate fea- 
ture has been established by extensive computer experiments and is illustrated here by the following 
interesting example. In Table 6b of Part I (Schwiderski, 1978a), one finds in the grid cell M = 236 
and N = 52 the obviously somewhat incoherent phase value 6 = 227°. This near-shore, deep-sea em- 
pirical datum (marked by a subbracket) was taken from publications by Irish et al. (1971), Now- 
roozi (1972), and others and was hydrodynamically interpolated into the tide model. The slightly 
anomalous datum was well noticed but retained as a possible retardation result of the nearby San 
Francisco Bay entrance and the proximity of the Farallon Islands. 


However, in 1979, Cartwright (see revised insert in Cartwright et al., 1979) communicated to 
the author the more plausible and correct phase value 6 = 206°, which he derived from the earlier 
original publication by Nowroozi et al. (1969). Evidently, the correct datum fits coherently into 
Table 6b of Part I (Schwiderski, 1978a,. where the erroneous value is interpolated) and into the new 
Tables (6M or 6N) in Appendix A, where the correct value is interpolated. 
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Naturally, it was mentioned in Schwiderski (1978a, 1979d, and 1980a) that the accuracy of 
the constructed tidal charts depends to some degree on the quality and quantity of the hydrody- 
namically interpolated data. After the arrival of the numerous and improved empirical data collec- 
tions mentioned above, the author found it justified to improve the constructed M, ocean tide by 
additional computations with revised empirical data inputs. 


In particular, the new data included slightly enhanced tidal amplitudes around the British Isles 
where earlier computations had to depend on marginal older data collections (see Schwiderski, 
1978a) because the tide tables of the British Admiralty (1977) listed no harmonic constants for the 
European waters. The new tape of tidal constants collected by the International Hydrographic 
Bureau (1978) filled this undesirable gap. Moreover, most of the numerous deep-sea tide data com- 
piled by Cartwright (see Cartwright et al., 1979) around the British Isles qualified as near-shore data 
and could be hydrodynamically interpolated into the new model (see Sect. 1, d. and h.). Similarly, 
the deep-sea tide data determined by the Canadian Hydrographic Service (see Cartwright et al., 
1979) for stations south of the Bay of Fundy between Nova Scotia and Cape Cod could be inter- 
polated as near-shore data, which led to slightly lower amplitudes. 


Finally, it must be mentioned that the phase data in the tide tables of the British Admiralty 
(1977) and on tape of the International Hydrographic Bureau (1978) differ by as much as 30° for 
the northwestern coast of India. Apparently, these unwanted differences can be traced to different 
interpretations of the time origin in both listings. After a telephone conversation with Mr. P. A. 
Bolduc of the Canadian Marine Environmental Data Service, who supervised the taping and revision 
of the data of the International Hydrographic Bureau, the data used earlier from the British Admir- 
alty were slightly revised. Subsequent computations indicated that the new phases seem to fit the 
general tide behavior in that area (e.g., Laccadiven tide data; Table 2M) somewhat better. 


A complete listing of all sources of empirical ocean tide data, which were interpolated into the 
new slightly improved M, tidal charts, is presented in Appendix A. In the meantime, Section 2 of 
this report lists the significant hydrodynamical input parameters that specified the reconstructed 
M, ocean tide. The major features of the global M2 tide are discussed in Section 3. A complete 
numerical display is presented in Appendix A where all tidal amplitudes and phases are gridwise 
tabulated in map-like charts. Corange (equi-amplitude) and cotidal (equi-phase) maps of the M, 
ocean tide are plotted in Appendix B. 


2. M; OCEAN-TIDE PARAMETERS 


The astronomical semidiumal principal lunar (M2 ) equilibrium tide 7 (or tide-generating poten- 
tial Gn; see Schwiderski, 1978a) at the geographical point (A, ¢) and instant (Y, D, t) is determined 
by 


n = K cos?@ cos(at + x + 2d) (1) 
where 

G = 9.81 m/sec? earth gravity acceleration 

dX = longitude (east in rad) 


@ = latitude (north in rad) 

Y 1975) = year number 

D = day number of year Y(D = | for January 1) 
t = universal standard time of day D (in sec) 


K = 0.242 334 m = M, equilibrium tide amplitude 
o = 1.405 19- 10~* sec”! = Mg tide frequency 
X = 2m(ho ~59)/180 = M, astronomical argument (in rad) 


= 279.696 68 + 36 000.768 930 4657+ 3.03 - 10°47? 

= mean longitude of the sun relative to Greenwich midnight of day D (in deg) 

= 270.434 358 + 481 267.883 141 377 -0.001 13377 + 1.9- 107° 7? 

= mean longitude of the moon relative to Greenwich midnight of day D (in deg) 
T = (27 392.500 528 + 1.000 000 035 6D] /36 525 

D = D+ 365(Y - 1975) + Int [(Y - 1975)/4] 

Int[x]) = integral part of x 
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The corresponding instantaneous ocean partial tide (Schwiderski, 1978a) is determined by 
¢ = &cos(or+x- 64), (2) 


where the local harmonic constants 


ore 
" 


&(A, ¢) = M2 ocean tide amplitude (in m) 


and 


mo 
W 


5, @) = M, ocean tide Greenwich phase (in rad) 
must be determined, say, by linear interpolation in the tidal charts of Appendix A. 


A simple second-order approximation in the sense of Love and Accad and Pekeris (see Part I. 
Schwiderski, 1978a, 1979c, and 1980; and Accad and Pekeris, 1978) yields 5 
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t€~0.612n and £20 = -0.0667¢, (3) 


ie., the corresponding terrestrial tide ¢° and the earth dip [© (yielding) under the oceanic tidal 
load {, respectively. A more elaborate and probably slightly more accurate earth dip £9 may be 
computed by using Farrell’s Green function (see Farrell, 1972 and 1973; and Schwiderski, 1980). 
In linear superposition, one finds the corresponding instantaneous geocentric partial M, tide: 


y= Fase 4 geo a) 


A detailed description of the hydrodynamical-empirical model to compute the ocean tidal am- 
plitudes & and phases 6 (listed in Appendix A) was given in Schwiderski (1978a, 1979c, d, and 
1980). In particular, all model input parameters such as the dimensionless eddy coefficient € (Eq’s. 
103 and 123), the bottom-friction paramater b (Eq’s. 4a and b), and the differencing parameters 
kK and « (Eq’s. 64 and 72) were all specified in Schwiderski (1978a) (referenced equations). These 
parameters were determined by extensive trial-and-error computations and are expected to remain 
unchanged for the construction of the leading 11 tidal components listed in Table 1 of Schwiderski 
(1978a). 


In the computation of the new M, tide model, the following mode-dependent parameters were 
used (see referenced equations in Schwiderski, 1978a): 


a. The time step At (Eq’s. 64, 123) 
At = 186.309 sec (5) 


b. The hydrodynamical interpolation control limits, k,, k,, and k3 (Eq’s. 88, 89, 94, 97, and 
99) 


k, = 0.045, k, =0.045, k3 =0.5. (6) 


It may be noted that the input parameters k, and k, of Equations 6 differ from k, = 0.06 and 
k, = 0.03, which were used (but erroneously listed in reverse order in Schwiderski, 1978a; see 
Schwiderski, 1979d) in the construction of the earlier M, tide model. The new parameter values 
were determined by additional trial-and-error computations to achieve best results. 


By comparing Tables 5-8 of Schwiderski (1978a) with the corresponding new Tables 8M and 
N, 6M and N, and 5M and N in Appendix A, one recognizes that the new parameters (which control 
the bottom-friction adjustments in the hydrodynamical interpolation) effected a change of about 
2 cm in amplitudes and about 3° in phases over some larger ocean areas. This slight change was con- 
sidered as a beneficial improvement relative to the improved empirical data. 


3. M2. OCEAN-TIDE FEATURES 


The entire constructed M, ocean tide is gridwise displayed in map-like amplitude and phase 
tables in Appendix A. The 41° x 71° charts cover the whole globe north of colatitude 169° (Antarc- 
tica) in three zones: a northern zone N from 0° to 71° colatitude, a middle zone M from 48° to 
118° colatitude, and a southern zone S from 98° to 168° colatitude. The overlapping geographical 
areas of the tidal charts have been chosen to provide a worldwide coverage for special applications 
and to allow the reader to scan the large amplitude and phase charts together in order to evaluate 
their quality and visualize the important tidal features. In addition, a generally superficial overview 
of some tidal features can be recognized by inspecting the more schematically plotted corange and 
cotidal maps provided in Appendix B. 


For an easy evaluation of the tidal charts in Appendix A, all hydrodynamically interpolated 
empirical tidal amplitudes and phases have been visibly marked by subbars for all shore data and 
subbrackets for all near-shore deep-sea input constants. Furthermore, the charts display the approxi- 
mate locations of distant off-shore deep-sea stations by subtildes under the computed amplitude 
and phase data. The corresponding empirical data, which were excluded from hydrodynamical inter- 
polation (see Sect. 1 and Schwiderski, 1978a, 1979d, and 1980), are listed and compared with the 
modeled data in Tables 1, 2, and 3. Finally, the approximate geographical locations of the impor- 
tant amphidromic points of zero amplitudes are marked by a circled (&). 


The tidal charts and maps permit the viewer to follow the tidal waves, that is the high water 
fronts (crests), in forward (or backward) direction, for instance, on their rotation around the amphi- 
dromic points. In the tidal phase charts of Appendix A, it is best to start from the prominently visi- 
ble 0° = 360° or 100° cotidal lines. Since the Greenwich phases specify the time lags (in degrees: 
30° = 1 hour) of the tidal crests relative to the cresting time of the corresponding equilibrium tide 
along Greenwich meridian (moon’s passage over Greenwich meridian; see Schwiderski, 1978a), one 
gathers a vivid impression of the significant global and local tidal phenomena. 


By following the tidal waves on their periodic rotations, one finds these waves passing through 
the specially marked stations in empirically correct time and with the correct height. In fact, all 
over the globe over 2 000 tidal phases and 2 000 amplitudes are coherently integrated. This is par- 
ticularly impressive for the charts of the Pacific Ocean, where the empirical data from so many clus- 
tered and scattered island stations fit smoothly into the surrounding computed tides. From the 
smoothness features of erratically interpolated tidal data (see Sect. 1), one concludes that this result 
is not an artifact of the interpolation applied but constitutes a vivid manifestation of the excellent 
compatibility of both the empirical and hydrodynamical procedures combined. 


On the basis of this observation, it can again (see Schwiderski, 1978a, c; 1979a, b, d; and 
1980) be estimated that the revised M, tidal charts permit a tide prediction with a uniform accur- 
acy of better than 5 cm anywhere in the open oceans. Naturally, near rough ocean basin reliefs (e.g., 
Arctic and Antarctic shores), where empirical tide (and depth) data are marginal in quality and 
quantity, a somewhat lesser accuracy must be expected. The estimated accuracy of the computed 


M, tide is, of course, fully validated by all 32 empirical tide data from distant off-shore deep-sea 
tide gauge stations, which are listed along with the computed data in Tables 1, 2, and 3. The differ- 
ences (not necessarily errors) range from 0 to 5 cm in amplitudes and 0° to 6° (12 minutes) in 
phases and thus verify the estimated prediction accuracy. 


From the tidal charts and maps in Appendixes 1 and 2, one concludes the following major fea- 
tures of the M, ocean tide. 


Table 1. North Atlantic Ocean Deep-Sea Empirical and Modeled M, Tides 


LATN | EMP—§ MOD€ JAPSONR SOURCES 


1.1.37 
1.1.29 
.1.30 
131 
1.32 
1.33 
1.34 
1.38 
1.39 
.1.40 
1.41 
1.42 
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75°38" 32°42’ 2. 3 
16° 25' 30°26’ 2 21) 
76°48" 28°27’ 8 2.15 
76°47 28°01’ 8 2.14 
67°32! 28°14’ 3 ee 
69°45’ 28°08’ 3 2. 4 
69°40’ 27°59’ 3 2. 8 
69°40’ 27°58’ 3 py ee 
69°20’ 26°28' 5 2.10 
69° 19’ 26°27' 5 2s 9 ’ 


€ = Amplitudes (cm) 

6 = Greenwich Phases (deg) 

IAPSO = Int. Assoc. for the Phys. Sci. of the Oceans 
C = Cartwright et al. (1979) 

M = Mofjeld (1975) 

P = Pearson (1975) 

Z = Zetler et al. (1975) 
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Table 2. Northeastern Pacific Ocean Deep-Sea Empirical and Modeled M, Tides 


LONG W EMP & EMP 65 MOD 6 IAPSONR SOURCES 


144°22' 56°08’ 
135°38' 53°19" 
132°47' 49°35' 
145°00' 34°00’ 
145°00' 34°00' 
124°26' 27°45’ 
129°01' 24°47’ 


MODE Aé 


& = Amplitudes (cm) 

5 = Greenwich Phases (deg) 

IAPSO = Int. Assoc. for the Phys. Sci. of the Oceans 
C = Cartwright et al. (1979) 

I = Irish et al. (1971) 

M = Munk et. al. (1970) 


Table 3. Southeast Indian Ocean Deep-Sea Empirical and Modeled M, Tides 


EMP 6 MOD6 _ Aé | IAPSONR__SOURCES 


132°01' 37°01' 4.1.1 C,IS 
132°09' 50°02' 4.1.2 C, IS 
132°07’ 60°01’ 4.1. 3 Cc, IS 


€ = Amplitudes (cm) 

5 = Greenwich Phases (deg) 

IAPSO = Int. Assoc. for the Phys. Sci. of the Oceans 
C = Cartwright et al. (1979) 

IS = Irish and Snodgrass (1972) 


a. Strong, Weak, and Incomplete Amphidromes. There exist amphidromic rotations around so- 
called ‘‘strong,”’ ‘“‘weak,”’ or “incomplete” amphidromes. The fully developed strong amphidromes 
display almost uniformly increasing amplitudes and phases in radial or tangential directions, respec- 
tively. An example of the strong amphidrome type may be seen in the North Atlantic (Table 9N). 
These amphidromes are distinctly isolated and independent of adjacent amphidromic rotations. 


Nr, sic 


Weak amphidromes are characterized by distinctly unevenly growing amplitudes and phases. 
These amphidromes occur mostly in pairs with a channel of low amplitudes and rapidly varying 
phases in between. A typical example can be seen in the northern South Atlantic (Table 9S). 


Finally, there exist areas with distinctly low amplitude tides and rapidly varying phases that 
seem to converge toward one (or more) points (incomplete amphidromes). Examples of this type 
may be seen in the North Pacific north of Midway (Table 5M) and in the center of the Caroline Is- 
lands (Table 4M). 


b. Compatible _and Incompatible Amphidrome Pairs. There exist pairs of amphidromes with 


“compatibly” (oppositely) or ‘“incompatibly” (parallely) rotating tidal waves. An example of the 
first kind, located in the South Atlantic, has a weak counterclockwise amphidrome in the center 
(Table 9S) and a clockwise amphidrome in the southeast (Table 1S). An example of the second kind 
exists in the North Atlantic. This pair has a strong counterclockwise amphidrome in the center and 
a neighboring one in the eastern Caribbean Sea. 


As a result of incompatible amphidromes, one finds lines in cotidal charts along which 
progressing tidal waves become standing and reverse their directions. On those wave reversal curves, 
cotidal lines are cuspate. Another interesting consequence of incompatible amphidromes is the pro- 
gressing tidal wave that follows along the whole eastern coast of the South and North Atlantic from 
South Africa in the south to Norway in the north. This tidal wave is obviously driven by the cone 
bined four amphidromic systems existing in the center of the South Atlantic (Table 9S), in the east- 
erm Caribbean Sea (Table 8M), in the central North Atlantic (Table 9N), and north of the Faeroe 
Islands (Table 9N), all of which rotate counterclockwise. 


c. Relief Distortions and _Retardations of Tides. As can be seen in all amplitude and phase 


charts of Appendix A, tidal waves crossing shallow continental shelves (see, e.g., the western Atlantic 
coast), narrow ocean ridges, or island chains (see, e.g., Caribbean, Hawaiian, Aleutian, or Marianas) 
suffer strong distortions and retardations. In fact, phase retardations of more than 30° and even 60° 
occur across ocean ridges. Moreover in certain coastal corners (e.g., in the North Sea, Table 1N; in 
the Bering Sea, Table 5N; on the Patagonian Shelf, Table 8S; and northwest of Australia, Table 4M), 
one finds oceanic tidal waves shedding small amphidromic rotations. As was pointed out in Schwi- 
derski (1978a, 1978d, and 1980), those rough tidal phenomena have long been known (see, e.g.. 
Harris, 1904; Defant, 1961; and Luther and Wunsh, 1975). The satisfactory resolution of those dis- 
tortions played a major role in the successful charting of ocean tides. 
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4. CONCLUSIONS 


The hydrodynamical interpolation technique has been applied to construct the semidiurnal prin- 
cipal lunar tide (M2) with an accuracy of better than 5 cm anywhere in the open oceans. The con- 
structed tide is displayed by tabulated charts in Appendix A and by corange and cotidal maps in 
Appendix B. The major features of the M, tide are discussed in Section 3. 


During the past three years, the author’s gradually improved M, tidal charts have been success- 
fully used in various applications. For example, Goad (1979) recently applied the M, tide model to 
reduce gravity measurements of the solid earth from oceanic tidal-load effects for stations in West- 
ern Europe, California, Japan, and Australia. In all cases, Goad estimated the accuracy of the com- 
puted gravity anomalies to be better than 4% microgals uniformly over all four areas. This important 
application provides another independent validation of the high accuracy of the utilized global M, | 
tide model. i | 
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APPENDIX A 


ATLAS OF 1” x 1° M, OCEAN TIDE 
AMPLITUDE AND PHASE CHARTS 
FOR 41° x 47° AREAS 


Bf ; APPENDIX A 
| ATLAS OF 1° x 1° M, OCEAN TIDE AMPLITUDE 
AND PHASE CHARTS FOR 41° x 47° AREAS 


1. GUIDE TO TIDAL CHARTS 


Amphidromic Points 
Subbars Mark Empirical Input Data at Shore Stations 
Subbrackets Mark Empirical Input Data at Near-Shore Deep-Sea Stations 


M = m: Longitude Number 

N = n: Colatitude Number 

Am = (m-0.5)°: Geographical Longitude East 
6, = (n-0.5)°: Geographical Colatitude 

Emjn = €Am, On): Amplitude (in cm) 

5m,n = 5(Am, 9p): Greenwich Phases (in deg. ; 30° = lh) 


Subtildes Mark Approximately Distant Offshore Deep-Sea Stations with Excluded 


Empirical Tide Data Listed in Appendix 2 


2. SOURCES OF EMPIRICAL TIDE DATA 
Publications: 


National Ocean Survey (1942), British Admiralty (1977), International Hydrographic Bureau 
(1978), Defant (1961), Miyazaki et al. (1967), Nowroozi et al. (1969), Munk et al. (1970), Zahel 
(1970), Irish et al. (1971), Irish and Snodgrass (1972), Nowroozi (1972), Luther and Wunsh (1975). 
Mofjeld (1975), Pearson (1975), Zetler et al. (1975), and Cartwright et al. (1979). 


Private Communications: 


D. C. Simpson (1977), National Ocean Survey, Rockville, Maryland; S. K. Gill and D. L. Porter 
(1978), National Ocean Survey. Rockville, Maryland; and K. Wyrtki (1978), University of Hawaii, 
Honolulu, Hawaii. 
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APPENDIX B 


ATLAS OF CORANGE AND COTIDAL MAPS 
OF THE M, OCEAN TIDE 


Amplitudes & of corange lines in 5, 10, 20, 30, 40, 50, 60, 80, 100, 120, 140, 160, 180, and 
200 cm. 


Greenwich phases 6 of cotidal lines in 15, 30, 45, 60, 75, 90, 105, 120, 135, 150, 165, 180, 
195, 210, 225, 240, 255, 270, 285, 300, 315, 330, 345, 360 = 0° where 30° = 1 hour. 
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Correction to Part II 
NSWC TR 79-414 


On page 5 correct h, to 


= 279.696 68 + 36 000.768 930 4857 + 3.03 + 10-4 7? 
= mean longitude of the sun relative to Greenwich midnight of day D (in deg) 


